Cubic Dresselhaus Spin-Orbit Coupling in 2D Electron Quantum Dots 
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We study eflects of the oft-neglected cubic Dresselhaus spin-orbit coupling (i.e., oc p^) in 
GaAs/AlGaAs quantum dots. Using a semiclassical billiard model, we estimate the magnitude 
of the spin-orbit induced avoided crossings in a closed quantum dot in a Zeeman field. Using these 
results, together with previous analyses based on random matrix theory, we calculate corresponding 
effects on the conductance through an open quantum dot. Combining our results with an exper- 
iment on conductance through an 8 ^rc? quantum dot [D. M. Zumbiihl et al, Phys. Rev. B 72, 
081305 (2005)] suggests that 1) the GaAs Dresselhaus coupling constant, 7, is approximately 9 
eVA'^, significantly less than the commonly cited value of 27.5 eVA'^and 2) the majority of the 
spin-fiip component of spin-orbit coupling can come from the cubic Dresselhaus term. 

PACS numbers: 73.21.La,73.23.-b,71.70.-d,05.45.Mt,73.63.Kv,72.20.My,72.20.-i 



Control over electron spin in semiconductors has 
promise for quantum computing and spintronics. In such 
applications, it is essential to understand how the trans- 
port of an electron through a circuit affects its spin; i.e., 
we must understand spin-orbit coupling. In technolog- 
ically important III/V semiconductor heterostructures, 
spin-orbit coupling originates in the asymmetry of the 
confining potential (called the Rashba term), which can 
be controlled by gates, and in bulk inversion asymme- 
try of the crystal lattice (called the Dresselhaus term). 
In quasi-2D systems, the Dresselhaus term has two com- 
ponents, one linear in the electron momentum and the 
other cubic. The cubic Dresselhaus term is usually ne- 
glected, as it is generally smaller than the linear contri- 
bution. Datta and Das proposed a spin-field-effect tran- 
sistor (SFET) for quasi-lD ballistic wires with Rashba 
coupling [l|. Schliemann, Egues, and Loss proposed an 
SFET that can operate in diffusive quasi-2D systems 
based on tuning the Rashba and linear Dresselhaus terms 
to be equal in strength, which produces long spin life- 
times, neglecting the cubic Dresselhaus term The 
strengths of the spin-orbit terms are difficult to measure 
independently, but a full understanding of their strengths 
is crucial to making such devices. Additionally, in con- 
fined systems such as quantum dots, some effects of the 
linear spin-orbit terms are suppressed Q, and it is im- 
portant to know the magnitude of the cubic Dresselhaus 
contribution, which could limit or even prevent the func- 
tioning of spintronic devices. 

We characterize the strength of the cubic Dresselhaus 
term in a confined system by its effect on avoided cross- 
ings in an in- plane magnetic field B|| that couples only 
to the electron spin. With no spin-orbit coupling, each 
eigenstate can be written as a product of orbital state \a) 
and spin quantized along By. Eigenstates \a t) and \(3 |) 
become degenerate when — e/5 = i?z, where ea,f3 are 
the orbital energies and Ez is the Zeeman energy pH . 
but spin-orbit coupling leads to avoided crossings. 

In the first half of this Letter, we estimate the cubic 



Dresselhaus contribution to the avoided crossings, which 
can be larger than the linear terms' contribution since the 
latter are suppressed for small B|| 0]. In the second half 
of this Letter, we relate these avoided crossings in closed 
quantum dots to the mean and variance of the conduc- 
tance when the quantum dot is connected to ideal leads. 
We compare these predictions to the results of Zumbiihl 
et al. on transport through an 8 /im^ quantum dot 0, Q 
and find that agreement is possible only if the cubic Dres- 
selhaus coupling constant in GaAs, 7, is considerably less 
than the frequently cited value of 27.5 eVA 0, Q from 
k • p theory. A smaller value of 7 has also been sug- 
gested by experiments [1, 0] and band structure studies 
[lol [III [1^ . Even with this smaller value of 7, we find 
that the cubic Dresselhaus term is the dominant spin- 
fiip mechanism in the sample considered. Both of these 
results are significant for spintronic devices. 

We consider conduction electrons in a 2D electron sys- 
tem grown on a (001) surface of a III/V semiconductor 
confined to a small area by a potential V{r). We use an 
effective Hamiltonian 
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where p = P — eA./ c, P is the canonical momentum, A 
is the vector potential from the perpendicular magnetic 
field, a is the vector of Pauli matrices, m is the effective 
mass. Ago = ei?ia-2/2Ai — e2?iCTi/2A2 is the effective spin- 
orbit vector potential, which contains both the Rashba 
and linear Dresselhaus spin-orbit terms, and Ai^2 are the 
(linear) spin-orbit lengths [l^, Q • We choose a coordi- 
nate system with axes ei=[110], e2=[liO], and e3=[001]. 
The second term is the cubic Dresselhaus term. 

In a system of linear size L, the linear spin-orbit terms 
can be gauged away to first order in L/\ by the unitary 
transformation Ti UHW = H' where U ~ exp(ir • 
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Aso) ^IS^I- Expanding to leading order in £/A, 

W = -^(p-ai-a||)2 + 6Z + 5Z 
Zm 

+ ^(p^-p?)(px^).e3 + nr), (2) 

where = -^^[e^i x r], a,, = -^^{xiGi/Xi + 
X20-2/A2)[e3 X r], 6^ = • (t/2, and b\ = 

-gUBiBixi/Xi + B2a;2/A2)tT3/4. 

When we apply a Zeeman field, we can treat each in- 
duced degeneracy as a two-level system, assuming the 
spin-orbit matrix elements, ego, are much less than the 
single spin mean level spacing, A — 2TTh^/mA, with m 
the conduction band efi'ective mass and A the dot area. 
The magnitude of the avoided crossings at the Fermi en- 
ergy is given by Cso = |(Q! Tl W° \P where ea-ep = Ez 
and Ect = Ef- We want to find the rms value of ego- Fol- 
lowing Ref. [sl, for a closed chaotic dot we may write 
A2 = (eso/A)2 as 
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FIG. 1: (color online) Normalized rms avoided crossing due 
to cubic Dresselhaus spin-orbit coupling, Ac, as a function of 
Zeeman energy Ez for four billiards with in-plane magnetic 

field along the indicated directions, with 7 — 8.5 eVA'^fl^; 

')aT:/3il ^i^a — — Ez)S{ea — Ep), (3) all billiards have A ^ 8 fim^ and n = 5.8 x 10^^ m"^. Insets 

show the billiard shapes with crystal axes. Solid lines indicate 
specular boundary conditions and dashed lines indicate diffuse 
boundary conditions, a) has a mixed phase space with small 
regions of regular trajectories, b) is a stadium billiard, c) is 
similar to the dot in Ref. [J , with diffuse boundaries to ensure 
chaos, and d) is a square with diffuse scattering from the top 
and bottom and specular scattering from the sides (see text). 



where the overbar indicates ensemble averaging and 
Oa,b^{a\0\b). 

We rewrite Eq. [3] as in Ref. [5] using the ^-dependent 
representation of the delta function and interaction pic- 
ture operators, and, after summing over f3, find 
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dt 



A27r?i 



a\rLf^{t)nf^\0)\a), (4) largest, highest density quantum dot in Ref. [J], with 



where loz — E^jfi is the Zeeman frequency, |a) is a 
typical orbital eigenstate with « iSp, and Tigp = 
li) is the spin-flip part of W° . We consider the 
cubic Dresselhaus term alone, = 2p'(P2 ~ V^){v\'^2 ~ 
P2fi), and estimate its contribution to A, which we call 
Ac. For example, if By = i3ei, 7i; 
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We estimate Ac semiclassically using a billiard model 
for the quantum dot, where the matrix element in Eq. [3] 
is replaced by the corresponding expectation value for a 
classical particle moving at the Fermi velocity, vp, start- 
ing at a random point moving in a random direction. We 
consider B± = for these simulations. Each of 2-3x10^ 
such trajectories is followed for an equal amount of time, 
which is generally about 300 bounces total in the for- 
ward and backward directions. Increasing the number 
of trajectories or bounces does not change the results. 
We calculate / dte-"^^^*-^">nyit)ny(t') for 100 ran- 
dom initial times t' on each trajectory as a function of 
ujz, and their average gives A^ when multiplied by the ap- 
propriate prefactors. We add a damping function to the 
integrand that sends it smoothly to zero as t approaches 
the simulation cutoff. 

We consider four different billiard shapes and, for 
specificity, choose parameters corresponding to the 



A = 8 iim^ and n = 5.8 x 10^'^ m"^. We use g = 0.44 
and m = 0.067too where rric is the electron mass. Fig. 
[Tfa-c) shows the resulting Ac{Ez) for these billiards, for 
three orientations of B||, with the Dresselhaus coupling 

constant chosen to be 7 = 8.5 eVA^. For other choices, 
Ac scales linearly with 7. The insets show the shapes of 
the billiards. For the semiclassical method of estimating 
Ac to be valid, we must have A <C Ez <C Et, where 
Et = tivp /VA is the Thouless energy. For the case dis- 
cussed here, A = 0.9 /ieV and Et ~ 80 fj,eV. 

We can understand the approximate scale of Ac by 
using a simpler, unphysical billiard. Consider an x Ly 
rectangle with specular reflection boundary conditions on 
the sides and diffuse boundary conditions on the top and 
bottom. At each collision with a diffuse wall, we choose 
the tangential momentum from a uniform distribution 
on [-pF, Pf] [IB- This choice gives a correct weighting 
for diffuse scattering and maintains detailed balance. In 
such a billiard, for B|| || x, in the limit Ez — > 0, 
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and we can break each trajectory into segments between 
collisions with the top/bottom walls. Along each seg- 
ment, p^ and Py are constant, and the particle takes time 
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t = mLy/ \py\ to move from one end of the segment to 
the other, so we can rewrite Eq. [5]as 
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which we can evaluate explicitly, since the j are uncor- 
related between segments. The particle begins moving in 
a random direction with V{px,o) = 7r~^(p|, — Px.o)~^^'^: 
where V is the probability density on [— pf, Pf]- Since 
the diffuse boundaries in this billiard are the top and bot- 
tom, T'{px.i^o) = ^/'2pF- We regularize the infinite sum 
by I]r=-oo(~l)" = 0, and, noting that Et = TipF/mLy, 
we find Ac = A'^'^p%, / {A^-n:'^ /^ETfi") . For the parameters 
in Fig. [H this gives IVc{Ez ^ 0) = 0.678. Finite values 
of Ez are not amenable to such simple treatment, but 
simulations of this billiard appear in Fig. [TJi, where the 
results for B|| || x, shown by the solid trace, approach 
the analytic prediction for Ez ~^ 0. 

Avoided crossings have not yet been directly measured 
in chaotic dots, but our calculations can be related to ex- 
periments measuring the conductance g through a quan- 
tum dot by Zumbiihl et al. 0, @]. To make this com- 
parison, we need a connection between avoided cross- 
ings in a closed dot and properties of the dot with leads 
attached. Cremers et al., using random matrix theory 
(RMT), worked out a similar connection for dots with 
only Rashba and linear Dresselhaus spin-orbit coupling 
|14j . We point out that the cubic Dresselhaus term can be 
added easily into the predictions of Cremers et al. without 
changing their formulas for {g) and var g by reinterpret- 
ing one of their RMT energy scales to include both linear 
and cubic spin-orbit terms. We now elaborate. 

In Ref. [14], the chaotic quantum dot is connected 
to two ideal leads with ^ 1 open channels, giving 
a scattering matrix from the circular orthogonal ensem- 
ble. They treat the magnetic field and spin-orbit coupling 
with a stub model [3] in which the stub has the M x M 
perturbation Hamiltonian, H'^y^-j,, given by association 
to H' in Eq. [2] (without the cubic Dresselhaus term), as 
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iAoixt + a±(T3) + za|| {Aicri + A2(J2) 



-b • (T + b^Bhcrs 



(7) 



where Ai, i — 0,1,2, are real antisymmetric matrices 
with {ti A^AJ) = S^j^P, Bh is a real symmetric matrix 
with (tr;B^) — M^, Al ^ 1 is the number of channels in 
the stub, and x, a±, a||, b, and b± are dimensionless pa- 
rameters, with X corresponding to B^, b to the Zeeman 
field, and a±, a||, and b± to the similarly named terms in 
Eq. [2] (without the cubic Dresselhaus term). Dephasing 
is included by setting Nctf = N + 2TTh/T^A, where is 
the dephasing time. Expressions are then obtained for 



(g) and var g as functions of x, aj 



b, b±, and N^ff 



to leading order in l/Nes [l3|- Zumbiihl et al. use these 
results to fit their data. 

Without the cubic Dresselhaus term, the correspon- 
dence between Eqs. [7] and [5] gives the following mapping 
from physical parameters to RMT parameters: 

7rKi;T/A($/$o)^ b^TTEz/A 
al = ttkEtIA{A/AXIo? (8) 
ajj - TiK'ET/A{AlA\lof{{L^l^if + {L2lMf) 
bl = ^n"{Ezf/ETA{A/A\lo) 

where $ is the magnetic fiux through the quantum dot, 
$0 — h/2e is the fiux quantum. Ago — %/ A1A2, Li^2 are 
the linear dimensions of the roughly rectangular dot, ori- 
ented along ei,2, a.nd k, k', and k" are geometric factors 
of order unity [l^, [l3|. We add the cubic Dresselhaus 
to this theory by noting that, as a random matrix, the 
cubic Dresselhaus term in Eq. [2] has the same symmetry 
as the a|| term in Eq. [71 i.e., it contains only ui and (72 
Pauli matrices. Assuming no correlation with the linear 
terms, we include the cubic Dresselhaus in Ti-'n^^-^ by set- 
ting a| = a| J -|- a| ^, where ayj is the Rashba and linear 
Dresselhaus contribution, given by Eq. [U and a||.c is the 
cubic Dresselhaus contribution. Since Ti^j^j^p contains the 
spin-orbit part of the Hamiltonian of the closed quantum 
dot, we relate ay^c to Ac by finding the rms spin-fiip ma- 
trix element (with spins quantized along By) due to a|[,cj 
giving a|| c = 27rAc. Including the cubic Dresselhaus term 
in TiJ^Mrp in this way cannot improve any fits to experi- 
mental data, as it does not change the fitting equations. 
However, its inclusion lifts the constraint that ay <C a 1 
[ist , similar to spatially varying spin-orbit strengths \u\ . 

Zumbiihl et al. observe weak anti-localization (WAL) 
in only one of the GaAs/AlGaAs heterostructure quan- 
tum dots they study [1, Q , and that dot gives the best 
defined values of the RMT parameters; we use it for the 
discussion of our results. The other dots do not contra- 
dict this discussion. The dot that displays WAL has area 
A — ^ fim^ and electron density n = 5.8 x 10^^ m~^, 
which are the values used in Fig. [TJ The dot has N — 2 
and A^cff= 13.9 [i]. 

Zumbiihl et al. measure var g as a function of By with 
time reversal symmetry broken by a small B±. They fit 
to the expression of Cremers et al. 1J|, with a|| (and all 
parameters except k") fixed to the value determined from 
the (g) data. We redo the fits to the var g data, constrain- 
ing only a|| > a||,c5 with a|| c from our simulations, and 
T0 fixed to the value determined from (g). From Fig.[Tl a 
typical value of Ac in all our billiard shapes is 0.4, giving 
a||,c = 2.5 (8.1) for 7 = 8.5 (27.5) eVA^ (recalling that 
Ac oc 7). We find that a value of oy c ~ 2.5 is compatible 
with the experimental data. However, if we require that 
a|| > 8.1, the fits to the data become markedly worse. 

Zumbiihl et al. also measure (g) as a function of B± , 
which they use to determine a||, finding a\\ = 3.1 0, [^, 
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Ml- Since of, 



if, „, we must have a\\ > a||^c: so we 



conclude that 7 = 27.5 eVA is inconsistent with these 
results, while 7 = 8.5 eVA is consistent with the data. 
So both (g) and var g data indicate that 7 should be 

closer to 9 eVA^than 28 eVA'^ (2^ . Moreover, even with 
the smaller value of 7, the cubic Dresselhaus term gives 
the dominant contribution to a\\. 

There are only a few other experiments pertaining to 
the value of 7 in GaAs. The best, most direct study is 
the Raman scattering in a GaAs/AlGaAs quantum well 
by Richards et al. in which they found 7 = 11.0 eVA Q; 
the same group also found 7 = 16.5 eVA in a different 
sample Q. A recent experimental value of 7 = 28 eVA 
from transport measurements [l9l | is less direct, includes 
the cubic Dresselhaus only as a density-dependent renor- 
malization of the linear Dresselhaus, and assumes the 
Rashba coupling is independent of gate voltage. Theoret- 
ical work has indicated that 7 is smaller in AlGaAs/GaAs 
heterostructures and superlattices than it is in bulk GaAs 
[lol . 11, [2^, so it is possible that experiments are not 
probing the bulk Dresselhaus coupling, though Ref. [l2| 
predicts 7 — 8.5 eVA in bulk GaAs. We include as sup- 



plementary information a table with experimental and 
theoretical values of 7 in GaAs and AlGaAs/GaAs. 

Strictly speaking, our calculations are not directly ap- 
plicable to the (g) data of Zumbiihl et al. , as our calcula- 
tions assume Ez S> A, and (g) is measured with B||=0. 
We do not believe, however, that ay changes signifi- 
cantly as Bii — !■ 0; similarly, Cremers et al. consider an to 



be constant for all Bm [14|. We believe that an f.(Bi| 



0) 



can be estimated by simply averaging our results from 
the different field directions in the limit By ^0. 



Our reinterpretation that = a| j + a| ^ requires, of 
course, that an j be less than 3.1 in the experiment of 



Zumbiihl et al. This reduction of 



can be absorbed 



into the geometric parameter k' (which was set to 1 with- 
out fitting in Zumbiihl et al.) without affecting any of the 
physical parameters, r^, Xso, found by Zumbiihl et al. 
Reducing k' is reasonable, as Ref. [l3| predicts k' ~ 1/3 
for a circular diffusive system. 

Since A cx A~^, and in a ballistic system Et oc 
^-i/2j^i/2^ we can see that if the thickness of the 2DEG 
does not change with density, a||j oc A'^^'^n^^^, while 
ay c (X Ac (X A^l^rv'^'^ . We therefore expect that the cu- 
bic Dresselhaus spin-orbit coupling should be relatively 
more important in small, high density dots, precisely the 
ones likely to be useful for producing an SFET. 

In summary, we have used billiard simulations to esti- 
mate the effect of the cubic Dresselhaus term on avoided 
crossings in a closed chaotic quantum dot. These results 
are related to the conductance through a dot with ideal 
leads attached. The cubic Dresselhaus plays a strong 
and previously ignored role in observed transport prop- 
erties through quantum dots. Our calculations suggest 



that 1) the Dresselhaus spin-orbit coupling constant, 7, 
in GaAs/AlGaAs heterostructures has a value near 9 
eVA and not the frequently cited value of 27.5 eVA , 
and 2) even with this smaller value of 7, in the exper- 
iments considered the cubic Dresselhaus term provided 
the bulk of the spin-fiip portion of the spin-orbit Hamil- 
tonian, which had previously been assigned to the effects 
of linear spin-orbit terms. The value of 7 in this techno- 
logically important system deserves further study. 
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the Fannie and John Hertz Foundation and NSF grants 
PHYOl-17795 and DMR05-41988. 
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TABLE I: Supplementary Information. Magnitude of Dres- 
selhaus spin-orbit coupling constant 7 in GaAs (eVA'^). All 
results are in bulk GaAs unless otherwise specified. 



Ref 



Technique /Comments 



EXPERIMENT 



[SI] 24.5 

[S2] 20.9 

[S3] 26.1 

[S4] 23.5 

[S5] 34.5 

[S6] 16.5±3 

[S7] 11.0 

[S8] 28±4 ; 31±3 



[S9] 



9 ; 8.5 



[SIO] 19 ; 30 
[Sll] 28 
[S12] 27.57 ; 21.29 



[S13] 14.0 
[S14] 14.9 ; 28.2 
24.12 
19.8 
8.9 



[S15] 
[S16] 
[S17] 



8.5 
24.21 



[S18] 
[S19] 

[S20] 10 ; 27.5 

[S21] 27.48 

[S22] 23.6 

[S23] 8.5 



Optical orientation using Dyalconov-Perel 
(DP) spin relaxation time. Conduction 
electrons in p-type GaAs. 
Same as [SI]. 

GaAs/AlGaAs heterostructure magneto- 
conductance measurements. Sets fc^ — > 
instead of (fc^). 

Raman scattering from I8OA GaAs quan- 
tum well (QW). Does not include Rashba 
term. 

Same as [S4]. No Rashba term. 

Raman scattering from GaAs/AlGaAs 

heterostructure. 

Raman scattering from asymmetric 
GaAs/AlGaAs QW. More data than 
earlier. 

GaAs/AlGaAs heterostructure magneto- 
conductance measurements. 

THEORY 

k • p theory ; Linear Muffin Tin Orbitals 
(LMTO). Room temperature band gap. 
8x8 ; 14x14 k ■ p models. 
14x14 k ■ p. 

k-p. Bulk GaAs ; Ga.65Al.35As/GaAs 

heterostructure. 

kp. 

LMTO ; k ■ p. 
5 level k ■ p. 
kp. 

spss* tight binding (TB) model of a 

GaAs/AlAs superlattice. 

TB model of 100 A GaAs/AlAs QW. 
5 level k ■ p. 
spss* TB model ; k • p. 
14 band 8x8 k ■ p theory. Notes that value 
is reduced to 19.6 in higher order pertur- 
bation theory (p. 74) 
TB model used to refine 14 band 8x8 k ■ p. 
Quasiparticle self-consistent GW method 
(QPscGW) 
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